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Abstract 


M molecular weight 


Effects of droplet interactions on drag, 
evaporation, and combustion of a planar droplet 
array, oriented perpendicular to the approaching 
flow, are studied numerically. The three- 
dimensional Navier-Stokes equations, with variable 
thermophysical properties, are solved using finite- 
difference techniques. Parameters investigated 
include the droplet spacing, droplet Reynolds 
number, approaching stream oxygen concentration, 
and fuel type. Results are obtained for the 
Reynolds number range of 5 to 100, droplet spacings 
from 2 to 24 diameters, oxygen concentrations of 
0.1 and 0.2, and methanol and n-butanol fuels. The 
calculations show that the gasification rates of 
interacting droplets decrease as the droplet spac- 
ings decrease. The reduction in gasification rates 
is significant only at small spacings and low 
Reynolds numbers. For the present array orienta- 
tion, the effects of interactions on the gasifica- 
tion rates diminish rapidly for Reynolds numbers 
greater than 10 and spacings greater than 6 droplet 
diameters. The effects of adjacent droplets on 
drag are shown to be small. 

Nomenclature 

B heat transfer number 


N$ number of species 

Nu Nusselt number 

p pressure 

Pr Prandtl number 

q iui ; heat flux 

r p droplet radius 
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ft universal gas constant 
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total drag coefficient 
friction drag coefficient, 
total friction force 
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specific heat; pressure drag coefficient, 

r total pressure force 
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x,y,z 
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Cartesian velocity component in x 
(axial )-direction 

Cartesian velocity component in 
y-direction 

Cartesian velocity component in 
z-di recti on 

mole fraction 

Cartesian coordinates 

mass fraction 


mass diffusivity of species i in the 
gas mixture 

droplet diameter 

total internal energy 

mixture fraction, defined as the mass 
originated from the droplets per unit 
mass of the gas mixture 


generalized curvilinear coordinates 

e oxygen concentration in the approaching 

stream; angle along the droplet surface 
measured from the front stagnation point 

u viscosity 

v stoichiometric coefficient (by mass) of 

oxygen 


enthalpy; heat transfer coefficient 


density 
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h° 

h f 

k 


heat of vaporization 
heat of formation 
thermal conductivity 


t stress tensor 

Subscripts 

A air property 


This paper is declared a work of the U.S. Government and is 
not subject to copyright protection in the United States. 


Copyright 1987 American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United Slates 
under Title 17, U.S. Code. The U.S. Government has a 
royalty-free license to exercise all rights under the copyright 
claimed herein for Governmental purposes. All other rights are 
reserved by the copyright owner. 


F 

fuel property 

m 

film condition 

prod 

combustion product property 

s 

droplet surface condition 

stoic 

stoichiometric condition 

UF 

unburned fuel property 

- 

approaching flow condition 

i 

index of species 


interest for practical sprays, by solving the 
three-dimensional Navier-Stokes equations for flows 
through droplet arrays. The arrays considered a^e 
monosized, planar, semi-infinite, with array-planes 
perpendicular to the approaching flow direction. 
(Another array orientation of practical importance, 
i.e., droplets arranged in tandem along the flow 
direction, is not considered in the present study.) 
To better simulate the flow around droplets, vari- 
able gas properties are used in the analysis. 
Numerical results are obtained for Reynolds number 
range of 5 to 100, droplet spacings of 2 to 24 
diameters, approaching flow oxygen concentrations 
of 0.1 and 0.2, and two types of fuel, methanol and 
n-butanol . 


Introduction 


Assumptions 


The evaporation and combustion of liquid fuel 
sprays have received considerable attention. Much 
of the theoretical work has focused on the trans- 
port of single, isolated droplets. j n regions 
near the fuel nozzle, however, dense spray effects 
are important and droplets may evaporate and burn 
quite differently from those simulated in the iso- 
lated droplet approach. 2 " 6 This is especially 
true for proposed advanced gas turbine combustion 
concepts, where liquid fuel and air are mixed in 
near stoichiometric proportions, in contrast to 
overall lean fuel/air ratios currently used in con- 
ventional combustors. The dense spray region is 
characterized 2 by atomization, droplet inter- 
action (defined here as modification of droplet 
transport rates due to the presence of adjacent 
droplets), droplet collision, coalescence, and 
breakup, change of turbulence properties by the 
droplets, the volume occupied by the liquid phase, 
etc. Droplet interactions are investigated numer- 
ically in the present study. 


Previous theoretical studies on droplet inter- 
actions have largely been limited to droplets in 
arrays or clouds in the absence of forced convec- 
tion (the diffusion theories), 6 *'^ although 
Stefan flow induced by evaporation may be included 
in the analysis. Calculations of this type indi- 
cated that interactions can significantly reduce 
droplet evaporation and burning rates even for very 
large droplet spacings. For example, for an 
array of four burning droplets at spacing of 10 
diameters, diffusion theory' predicted that the 
droplet lifetime increased by 20 percent over the 
lifetime of a single, isolated droplet. For the 
same droplet spacing, in the presence of forced 
convection, the effects of interactions would be 
negligible, even for a droplet Reynolds number as 
low as two. 9 ’* 0 Most investigators also adopted 
the constant property assumptions which resulted in 
predicted flame sizes exceeding the experimental 
values by factors of three to five.* Because 
larger flames compete more extensively for oxygen, 
the constant-property models predict much stronger 
interactions compared to the experimental data.^ 
Since most practical sprays involve appreciable 
droplet Reynolds numbers and large variations of 
thermodynamic and transport properties in the flow 
field, the constant-property diffusion theories 
appear to have limited utility in the analysis of 
droplet interactions in combusting fuel sprays. 


The present study considers droplet inter- 
actions in a steady-state situation in the presence 
of forced convection, covering Reynolds numbers of 


The equations describing the flow include the 
conservation of mass, momentum, energy, and species 
for both the gas and the liquid phases. Additional 
constitutive relations are the equation of state 
and the thermodynamic and transport properties as 
functions of temperature, pressure, and species 
concentrations. In order to make the problem 
tractable and to avoid undue complications, the 
following assumptions are made: 

(1) The gas-phase processes are quasi-steady, 
i.e., the qas phase adjusts to the steady state 
structure for the imposed boundary conditions at 
each instant of time. This assumption is justified 
by the large liquid/gas density ratio. Because of 
the density difference, the liquid-phase proper- 
ties, e.g., the surface regression rate, surface 
temperature, and species concentrations, change at 
rates much slower than those of the gas-phase 
processes. 

(2) Liquid-phase internal motion and transient 
heating are neglected. The droplet remains uni- 
formly at the wet-bulb temperature, which is deter- 
mined by balancing the heat transfer to the liquid 
and the latent heat of vaporization. For single- 
component fuel at low or moderate ambient pres- 
sures, Law and Sirignano^ have indicated that 
transient heating constitutes only a small fraction 
of total interphase energy transport after the 
initial 10 to 20 percent of the droplet lifetime. 

Because of assumptions (1) and ( 2 ), the drop- 
let size, spacing, and Reynolds number all remain 
at their initial values, and the transient heating 
effect is excluded from the analysis. The assump- 
tions offer two very important advantages. First, 
the detailed flow field solutions within the liquid 
droplet become unnecessary, and the only interphase 
properties sought in the solution procedures are 
temperature (wet-bulb temperature) and pressure, 
from which other interphase boundary conditions 
(such as fuel vapor concentration) can be calcu- 
lated. Second and more important, the parametric 
effect of individual factors affecting droplet 
interactions can then be studied, without complica- 
tions caused by simultaneous change of more than 
one parameter and the history of evaporation. 

(3) Phase equilibrium is maintained at the 
droplet surface. The fuel vapor concentration at 
surface is given by the saturated vapor pressure 
correlation for the oure liquid, e.g., the Clausius- 
Clapeyron equation, at the wet-bulb temperature. 
Surface tension corrections are neglected. 
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(4) Effects of thermal radiation, turbulence, 
and the Dufour and Soret effects are neglected. 

The droplets are spherical in shape. The pressure 
of the approaching flow is maintained at 1 atm and 
the ambient gases have negligible solubility in the 
liquid phase. The effect of natural convection is 
neglected since the Grashof number is generally two 
orders of magnitude smaller than the Reynolds 
number for the flows considered in this study. 

(5) Mass diffusion is represented by an effec- 
tive binary diffusion law. 

(6) The chemical reaction rates are much 
faster than the gas-phase transport rates such that 
combustion occurs at a thin flame sheet where fuel 
and oxygen meet in stoichiometric proportions. The 
combustion processes proceed to completion. 
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( 1 ) 


(7) The mass diff usi vities of the combustion 
products and the fuel vapor are equal. Because of 
assumptions (6) and (7), the gas-phase species con- 
centration field of the burning droplets can be 
described by one conserved scalar quantity, e.g., 
the mixture fraction (defined as the fraction of 
mass originated from the droplets). 

(8) Effects of the wake instability are 
neglected. The onset of wake instability for solid 
particles in isothermal flow occurs at Reynolds 
number around 130, which is greater than the maxi- 
mum Reynolds number (100) considered in the present 
study. No information is available concerning 
Reynolds number for the onset of wake instability 
for vaporizing or burning droplets. 
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and U = £ x u + E y v + e z w , 

V = n x u + n y v + n z w , 
W = Cx u + 5y v + Cz w 


( 2 ) 


(3) 


The above assumptions may become invalid at 
very low pressures (much below 1 atm), near thermo- 
dynamic critical point of the fuel, for very small 
droplets (on the order of 1 ym), or in the presence 
of luminous flames. 

Analysis 

Governing Equations 

The three-dimensional, unsteady Navier-Stokes 
equations are solved for the asymptotic steady- 
state flow field in droplet arrays. The equations 
are cast in conservation law form and solved using 
a finite-difference method. To enhance numerical 
accuracy and efficiency, coordinate mappings are 
employed which bring droplet surface and symmetric 
planes onto coordinate surfaces, and cluster grid 
points near the droplet surface. The governing 
equations, written in the generalized curvilinear 
coordinates c(x,y,z). n(x,y,z), and c(x,y,z), 
are given as follows: 12 


where U, V, and W are contravariant velocities, 
S x , Ky, K Z y et c. are the metric coefficients and J 
is the Jacobian of the coordinate transformation. 

The viscous flux terms are given by 
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* The forms for F v and G v are similar to 
E v , except K is replaced by n and c in F v 
and Gi v , respectively. 
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The stress and viscous dissipation terms and 
the species and thermal energy diffusion terms are 
gi ven by 
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where the subscripts x, y, and z denote differ- 
entiation in the respective directions. The total 
internal energy and pressure are given by 


e = X! ph - Y . - p + £(u 2 + v 2 + w 2 ) (6) 

i=l 1 1 l 



. s Y- 

and P = pftT X ^ (8) 

where h°. is the heat of formation for species 

i at the reference temperature Tr. The 
Cartesian derivatives are to be evaluated in 
space via the chain-rule, for example 


T x “ ^x T c + n x T n + ^x T c 
Thermodynamic and Transport Properties 


The specific heat, thermal conductivity, and 
viscosity for each species are determined by poly- 
nomials of temperature, such as, 

C Pi = A. + B.T + CT 2 + D^ 3 (9) 


The coefficients of these polynomials are found in 
Ref. 13. The specific heat of the gas mixture is 
obtained by concentration weighting of each 
species. The thermal conductivity and viscosity 
of the mixture, however, are calculated using 
Wilke's law, 14 for example, the mixture viscosity 
is determined by 
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where 



(U) 


The binary mass diffusivity for the fuel vapor 
in the ambient gas is obtained using the Champman- 
Enskog theory in conjunction with the Lennard-Jones 
i ntermolecul ar potential-energy functions. Details 
of this method can be found in Ref. 14. 


Combustion Model 


Both droplet evaporation and combustion are 
considered in the present study. For the burning 
droplet case, a mixing controlled combustion model 
is employed. Chemical reaction rates are assumed 
to be much faster than the gas-phase mixing rates, 
and the chemical reactions proceed immediately to 
completion when the fuel vapor and the oxidizer are 
mixed in stoichiometric proportions. If we further 
assume that the combustion products have the same 
binary mass diffusivity as the fuel vapor in the 
gas mixture, the flame front positions can be deter- 
mined from the stoichiometric mixture fraction 
values. The concentrations of the unburned fuel 
vapor, combustion products, oxygen, and nitrogen 
(assuming the approaching flow is composed of only 
oxygen and nitrogen) can then be determined by the 
mixture fraction and the flame front location. 
Denoting the stoichiometric coefficient (by mass) 
of oxygen as v and the oxygen mass concentration 
in the approaching flow as e, the stoichiometric 
mixture fraction value and the species concentra- 
tions in the gas mixture can be calculated from 

f = - 

stoic i+l 

1 e 

from the droplet surface to the flame front 

H) 

[Y n2 ] = (1 - e) (1 - f) 

[y 02 ] = °- 

[Y up ] = 1 - [Y prod ] - [Y n2 1 - [Y 02 ] 


( 12 ) 
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and from the flame front to the outer edge of the 
flow domain 

C Y prod ] " f(1 + V) 




(1 

0. 


[ y 0 2 ] = 1 


e) (1 - f) 


tY prod^ tY N2^ 




(14) 


After the concentration field is obtained, temper- 
atures and pressures are calculated from Eqs. (6) 
to (8), using Newton's iteration method. 

Surface Integral Parameters 

Previous numerical and experimental studies 
on drag and heat and mass transport for isolated 
droplets in high-temperature flows are abundant .I 5 
Results from these studies are used to validate the 
analysis and the numerical method described in this 
paper. Since results were presented in the form of 
drag coefficients and Nusselt numbers for most of 
the existing studies, these integral parameters are 
also calculated in the present study to facilitate 
comparison. 

The drag force on the liquid droplet consists 
of contributions from the viscous stresses, the 
pressure, and the momentum flux at the interface. 
The computed momentum flux force (the thrust drag) 
at the droplet surface is about two orders of mag- 
nitude smaller than the other two forces and is 
therefore neglected. If grid orthogonality is 
maintained at the surface, the axial (approaching 
flow direction) component of the surface shear 
stresses, can be written, in terms of variables in 
the curvilinear coordinates (s,n,t), as 
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the rate ot dilation. 

The expressions for g nn and g cc are similar 

to g^. 

The axial component of the pressure force is 
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Integrating over the droplet surface and 
nondimensional izing with approaching flow 
quantities, the drag coefficient becomes 
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To be consistent with most of the published data, 
the Nusselt and Praridtl numbers are calculated 
using the film properties, i.e.. 
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where the subscript m refers to the film condition 
defined by e = 1/2 in the following equations 

T m = 0 T S + (1 - e) ^ for evaporating droplets 

Y m ® 0 T S + (1 - ©) Aflame for burning droplets 

Y m = eY s + (1 - e) Y ® for evaporating droplets 

and 

Y m = eY s + (1 - e) Y flame for burning droplets (18) 


The heat transfer coefficient is given as 

for evaporating droplets 
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for burning droplets 
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where the derivative a/ar, for orthogonal grids 
at the surface, is given by 
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Numerical Methods 

Grid System. Semi-infinite planar arrays of 
equally spaced droplets are employed in the present 
study. A schematic of a typical array configura- 
tion is shown in Fig. 1. Because of the symmetric 
arrangement of droplets, only a 45° sector, as 
indicated in Fig. 1, needs to be considered in the 
computation. To enhance numerical accuracy and 
efficiency, coordinate mappings are used which 
bring droplet surface and symmetry planes onto 
coordinate surfaces, and cluster grid points near 
the droplet surface. This would also help the 
implementation of boundary conditions, since no 
interpolations are required at boundary surfaces. 

An 0-type grid, as shown in Fig. 2, is generated 
algebraically, with minimum radial spacing (in the 
physical domain) of 0.02 droplet radius, and the 
grids are stretched exponent! al ly in the radial 
direction outward from the droplet surface. For 
clarity of presentation, much larger grid spacings 
near the droplet surface and fewer grid lines are 
shown in Fig. 2 than actually used in the calcula- 
tions. The grid orthogonal i ty is maintained at and 
near droplet surface so that the interphase heat 
and mass fluxes and the shear stresses can be more 
easily calculated. In the curvilinear coordinates 
the computational domain is rectangular parallel- 
piped with uniform grid spacing, which facilitates 
the use of standard unweighted differencing schemes 
and helps to maintain higher order numerical accu- 
racy. Most of the calculations in the present 
study are performed using a 55 by 15 by 55 grid, in 
the axial (the approaching flow direction), azimu- 
thal, and radial directions, respectively. Fifteen 
grid lines in the azimuthal direction are consid- 
ered adequate for spatial resolution, since the 
flow domain (only 45 ) and the gradients of flow 
properties are smaller in this direction than in 
the other two directions. The grid is reduced to 
41 by 12 by 41 for the lowest Reynolds number case 
(Re = 5) due to the numerical instability in the 
finer grid. One calculation for Re = 100 using a 
65 by 18 by 65 grid is carried out and the solu- 
tions show negligible improvement over the results 
obtained using the 55 by 15 by 55 grid. Therefore, 
it may be concluded that the numerical solutions 
are relatively independent of the grid density. 

Finite-Difference Procedure 


although much easier than the unfactorized opera- 
tors, accounts for the major portion of the total 
computational effort of the implicit scheme. To 
improve the numerical efficiency a similarity 
transformation 19 is employed, which diagonalizes 
the blocks in the implicit scheme and produces 
scalar pentadiagonal operators in place of the 
block operators. 

Boundary Conditions 

The boundary conditions are implemented expli- 
citly. The velocity, temperature, static pressure, 
and species concentrations are specified for the 
approaching flow. At the downstream plane where 
the flow leaves the computational domain, flow 
properties are extrapolated from interior points 
except for the static pressure, which is set equal 
to the approaching flow value. These upstream and 
downstream boundary conditions are applied at a 
distance of 25 diameters from the center of the 
droplet. At the mid-planes between the droplets, 
symmetry conditions are applied. The droplet sur- 
face mass flux due to gasification is given by 



and the gas velocity components at the surface are 
obtained accordingly. The pressure on the droplet 
surface is calculated with a normal momentum rela- 
tion (obtained by combining the three momentum 
equations) A* The droplet surface temperature 
is taken as the wet-bulb temperature, which is 
obtained from the balance of the total heat trans- 
fer to the surface and the latent heat of vapori- 
zation and, therefore, is part of the solution. 

The surface fuel vapor concentration Yp s is 
obtained from the partial pressure of the saturated 
fuel vapor at the wet-bulb temperature, using the 
Clausius-Cl apeyron equation. For the case of drop- 
let evaporation (nonburning), the gas phase is con- 
sidered as a binary mixture of fuel vapor and air, 
and the mixture fraction f is equivalent to the 
fuel vapor concentration; hence, f s = Yp s and 
Y a$ = 1 - Yp s at droplet surface. For burning 
droplets, the surface mixture fraction is 
calculated by 


The finite-difference scheme used for solving 
the governing equations is the delta-form implicit 
approximate factorization algorithm described by 
Beam and Warming. 16 Since this and other 
similar schemes are fully documented in the 
literature, 1 ^* 16 * 1 ' only a very brief discussion 
of the numerical method will be given here. 

Because only the asymptotic steady state solutions 
are required, a first-order Euler implicit scheme 
is used to integrate the unsteady Navier-Stokes 
equations in time. The spatial derivative terms 
are approximated with fourth-order central differ- 
ences. Fourth-order explicit and second-order 
implicit artificial dissipation terms are added to 
the basic central -differencing algorithm to control 
the nonlinear numerical instabil ity. 16 Local time 
linearizations are applied to the nonlinear terms 16 
and an approximate factorization of the three- 
dimensional implicit operator is used to produce 
locally one-dimensional finite-difference oper- 
ators. 16 * 17 The resulting operators are block 
pentadiagonal matrices, and their inversion. 


f s 


+ f stoic^ 


( 21 ) 


where f s toic is 9 iven b Y Ed* ( 12 ) . Concentra- 
tions of the remaining species at the surface are 
obtained using Eq. (13) . 


Results and Discussion 


Calculations are first made for single, iso- 
lated solid particles and evaporating droplets, 
where the abundance of existing experimental and 
numerical data facilitates the validation of the 
analysis described in previous sections. The 
streamlines and isotherms for an isolated methanol 
droplet in a hot air stream at a temperature of 
800 K, a Reynolds number of 100, and a pressure of 
1 atm are shown in Fig. 3. There is a recircula- 
tion region (wake) formed behind the droplet and 
the reattachment point at the axis (e = it) is 
located at 0.96 diameters from the rear stagnation 
point. The isotherms show steep gradients near the 
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front stagnation point, indicating that the heat 
transfer rate is higher at the front half of the 
sphere. The locations of the reattachment points 
behind solid particles in isothermal flows are 
also calculated for Reynolds numbers ranging from 
20 to 100. The predictions (not shown here) agree 
very well with the calculations by Rimon and 
Cheng^O and the experimental data quoted by the 
same authors. As an example, the present calcula- 
tions show that the reattachment distance for a 
solid particle at Re of 100 was 0.92 diameter, 
compared to the value of about 0.90 diameter 
reported in Ref. 20. 

The drag of isolated solid particles in iso- 
thermal flows, and the drag and heat transfer of 
isolated evaporating droplets in hot streams at a 
temperature of 1000 K are compared with the numer- 
ical results reported by Renksizbulut and Yuen. 

The friction, pressure, and total drag coefficients 
are shown in Fig. 4, and the heat transfer results 
are shown in Fig. 5. The agreement between the 
present calculations and the results of Ref. 15 is 
very good. Since the numerical results in Ref. 15 
correlate well with a wide range of experimental 
data, the present numerical results also are in 
good agreement with experimental data. Figures 4 
and 5 also indicate that the standard drag law and 
the conventional empirical expressions for Nusselt 
number can be used for evaporating droplets in 
flows with large variations of transport proper- 
ties, provided the proper film properties are 
chosen for evaluation of Reynolds number and the 
heat transfer number (B). 

The favorable comparisons discussed above have 
demonstrated the validity of the present analysis 
and numerical procedures. Therefore, we can pro- 
ceed with confidence with the calculations of the 
interacting droplets. The droplet assemblages 
considered are planar arrays of equally spaced 
monosized droplets. The arrays are oriented per- 
pendicular to the approaching flow direction. In 
the following, the effects of interactions are 
presented as a gasification rate correction factor 

_ rate of gasification of a droplet in an array 
" rate of gasification of an isolated droplet 

The gasification rate correction factors for 
evaporating (nonburning) methanol droplets are 
shown in Figs. 6(a) and (b) for T^ = 700 K 
and 1400 K, respectively. As seen in Fig. 6, 
droplet interactions are only important for small 
spacings and low Reynolds numbers. They become 
negligible for spacings greater than about 6 dia- 
meters and Reynolds numbers greater than about 10. 
The present calculations show much weaker and 
shorter-ranged interactions than predicted by the 
diffusion theories, 7 ’ 8 where the effect of forced 
convection is not considered. A close inspection 
of the predicted flow field indicates that, in the 
presence of forced convection, temperature and 
concentration variations are contained in a thin 
boundary layer around the droplet, and the 
approaching stream conditions prevail outside this 
boundary layer. Since the boundary layer thickness 
around the droplet is of the order of magnitude of 
one droplet diameter for the Reynolds numbers con- 
sidered here, the effects of neighboring droplets 
on evaporation are not likely to be very signifi- 
cant for droplet spacings much greater than one 
diameter. 


The results of Figs. 6(a) and (b) are very 
similar, except Fig. 6(b) shows slightly stronger 
interactions. The thicker thermal and concentra- 
tion boundary layers and the stronger competition 
among neighboring droplets for thermal energy 
caused by the more intense evaporation at higher 
temperature are responsible for the increased 
interactions. 

Law et al.l> 6 have indicated that the flame 
size and the ambient oxygen concentration are the 
major factors that determine the extent of inter- 
actions of burning droplets. The flame shapes for 
single, isolated methanol and n-butanol droplets 
at Re = 25 are illustrated in Figs. 7(a) and (b), 
respectively. Two approaching flow oxygen concen- 
trations, i.e., 0.1 and 0.2, are considered. Since 
the stoichiometric mixture fraction of n-butanol 
fuel is smaller than that of the methanol fuel, the 
flame size of the former is larger than the latter, 
especially near the wake where the high mixture 
fraction region extends to several diameters down- 
stream of the rear stagnation point. The influence 
of oxygen concentration on flame size is also clear 
from the figures. The stoichiometric condition 
occurs farther away from the droplet surface for 
the lower oxygen concentration flow, yielding 
larger flame stand-off distance. 

Another important aspect in droplet transport 
analysis is the evaluation of physical properties 
for the gas mixture around the droplets >5,15 
To illustrate the importance of physical proper- 
ties, the gas-phase thermal conductivity, viscos- 
ity, and the product of density and fuel vapor mass 
diffusivity in the e = 90° plane, normalized by 
the interphase properties, are plotted against 
radial distance in Figs. 8(a) and (b) for an iso- 
lated droplet undergoing evaporation (Fig. 8(a)) 
or burning (Fig. 8(b)). The figures clearly show 
that, if physical properties are taken to be con- 
stant and evaluated at free stream (evaporation 
case) or flame (burning case) conditions, droplet 
transport rates will be significantly overestimated 
compared to the variable-property approach. Law 
et al. 1 * 6 have also pointed out that, for iso- 
lated droplets, because of constant property 
assumptions, theoretical predictions of the flame 
size consistently exceed the corresponding experi- 
mental values by factors of three to five. 

The gasification rate correction factors for 
burning methanol droplets are shown in Figs. 9(a) 
and (b), for approaching flow oxygen concentrations 
of 0.1 and 0.2, respectively. The interactions are 
stronger than the results shown in Figs. 6(a) 
and (b) for the evaporating droplets. This is 
attributed to the higher gasification rate due to 
the presence of the flame and the competition for 
oxygen by the neighboring flames in the burning 
case. Figure 9(a) also shows stronger interactions 
than those in Fig. 9(b), which can be explained on 
the basis of flame size in that the larger flame in 
lower oxygen concentration stream competes more 
vigorously for oxygen with neighboring flames and 
tends to have stronger interactions. 6 It can be 
clearly seen in both Figs. 9(a) and (b) that, in 
contrast to the findings based* on the diffusion 
theories, 7 * 8 the effects of interactions diminish 
rapidly for droplet spacings greater than 6 dia- 
meters and Reynolds numbers greater than 10. The 
stronger and longer ranged interactions obtained 


7 


by the diffusion theory are mainly due to the 
absence of the forced convection and the use of 
constant properties in the analysis. 

The burning rate correction factors for 
n-butanol droplets are shown in Figs. 10(a) 
and (b). The differences between Figs. 10(a) 
and (b) are similar to those between Figs. 9(a) 
and (b), i.e., larger flame size in lower ambient 
oxygen case yields stronger interactions. 

Figures 10(a) and (b) show that, compared to 
Figs. 9(a) and (b), n-butanol droplets experience 
slightly stronger interaction effects than methanol 
droplets under same flow conditions. This can 
again be explained on the basis of flame size in 
that the larger flames of the n-butanol droplets 
(as shown in Figs. 7(a) and (b)) compete more vig- 
orously for oxygen since they are physically closer 
than the smaller flames of the methanol droplets. 

The effects of interactions on droplet drag 
are also investigated in the present study. The 
blockage of flow by the adjacent droplets acceler- 
ates the flow (venturi effect) and produces larger 
shear stresses as well as a larger wake, and 
thereby a slight increase in friction drag and 
pressure drag. The increase in drag due to the 
venturi effect, however, is somewhat mitigated by 
the reduction in the boundary layer viscosity due 
to the lower flow temperatures around the inter- 
acting droplets (resulted from the competition for 
thermal energy for the evaporating droplets and the 
larger flame stand-off distance for the burning 
droplets). The net result is an insignificant 
change of drag due to interactions. 

Conclusions 

In the present study, we have investigated the 
effects of droplet interactions on the drag and 
gasification rates of evaporating and burning drop- 
let arrays. The arrays considered are planar, 
semi-infinite, and are composed of equal-spaced 
monosized droplets of same fuel type, with array 
planes perpendicul ar to the approaching flow direc- 
tion. The following conclusions can be drawn. 

1. The present analysis predicts less intense 
and much shorter-ranged interaction effects than 
those obtained using the diffusion theories. The 
difficulties with the diffusion theories lie in the 
fact that forced-convection and variable-property 
effects are neglected in the analysis. Since most 
practical sprays have appreciable droplet Reynolds 
numbers and involve large property variations, the 
present analysis appears to have more relevance in 
the consideration of droplet interactions compared 
to the diffusion theories. 

2. The effects of droplet interactions are 
stronger for the type of fuels and ambient oxygen 
concentrations which allow for larger flame sizes 
and thereby more intense oxygen competition. In 
this case, low ambient oxygen concentration and the 
fuel with smaller stoichiometric fuel mass fraction 
are in favor of interactions. 

3. The effects of droplet interactions on 
drag are small for the array configuration con- 
sidered in the present study. 

Arrays with droplets aligned in tandem along 
the flow direction are not considered in the 


present study. Since, for droplet Reynolds numbers 
encountered in practical sprays, the wake can 
extend to more than one droplet diameter downstream 
of the rear stagnation point and the wake flame is 
even longer, the interaction effects for tandem 
arrays are likely to be significant and warrant 
further study. 
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(a) Vertical plane (x = 0). 



(b) Meridional plane (y = 0). 
Figure 2. - Grid system. 





Figure 3. - Streamlines and isotherms of a vaporizing methanol droplet 
for Too * 800 K and Re * 100 at meridional plane (y - 0). 
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(r - r p )/r p 


(b) Burning droplet, Yg^ = 0-2, Too = 400 K - 

Figure 8. - Normalized gas mixture transport properties for 
isolated methanol droplet at 0 = 90° plane, Re » 50. 
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